Anomalous particle-number fluctuations in a three-dimensional interacting 

Bose-Einstein condensate 
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The particle-number fluctuations originated from collective excitations are investigated for a three- 
dimensional, repulsively interacting Bose-Einstein condensate (BEC) confined in a harmonic trap. 
The contribution due to the quantum depletion of the condensate is calculated and the explicit 
expression of the coefficient in the formulas denoting the particle-number fluctuations is given. The 
results show that the particle- number fluctuations of the condensate follow the law < SNq >~ 
jy22/i5 ^ e fluctuations vanish when temperature approaches to the BEC critical temperature. 
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The remarkable experimental realization of Bose- 
Einstein condensation in trapped, weakly interacting 
atomic gases has stimulated intensive theoretical and ex- 
perimental researches on Bose condensed gases 0. In 
particular, much attention has been paid to the study 
of the particle-number fluctuations of Bose-Einstein con- 
densates (BECs)0. This is mainly due to the fact that 
the particle-number fluctuations play a central role in 
the understanding of the statistical properties of BEC. 
In addition, the particle-number fluctuations of a con- 
densate may result in the fluctuations of chemical poten- 
tial and thus lead to the phase diffusion of condensate Q. 
Therefore, the study on the particle-number fluctuations 
of BEC is not only of an intrinsic theoretical interest 
but also useful to understand and control the coherent 
property of BEC. On the other hand, a clear theoretical 
understanding will be helpful to guide a direct experi- 
mental observation on the particle-number fluctuations 
in a BEC. 

Up to now there exist a lot of theoretical works ex- 
ploring the property of the particle-number fluctuations 
in BECs. In an ideal Bose gas the particle- number fluctu- 
ations < SNq > =< Nq > - < N > 2 have been studied 
rather thoroughly in a homogeneous case (e. g. in a 
box)0 and an inhomogeneous case (e. g. in a trap)0. 
The role of interatomic interaction in the particle-number 
fluctuations is an important theoretical problem and 
hence intensive theoretical investigations have been given 
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the particle-number fluctuations is described by the value 
of the power 7 in the expression < SNq >~ iV 7 with N 
being the total number of particles. Refs.0,0,01 give 
the result of anomalous fluctuations with 7 = 4/3, while 
Ref. argue that the fluctuations should be normal with 
7 = 1. 

For the temperature far below the critical tempera- 
ture of BEC, the collective excitations created from con- 
densate play an important role in the particle-number 



fluctuations. The physical reason is that due to the in- 
teratomic interaction the creation of the collective ex- 
citations induces a change of the particle number in 
the condensate. In the last few years, the particle- 
number fluctuations originated from collective excita- 
tions have been studied for three-dimensional (3D) 0,0, 
and 2D 0| weakly interacting Bose-condensed 
gases. In a pioneering work by Giorgini et al 0, 
the particle-number fluctuations are investigated within 
a traditional particle-number-nonconserving Bogoliubov 
method. Kocharovsky et al0] extended the results of 
by using a particle-number-conserving operator formal- 
ism. The scale behavior of the interacting condensate in 
a box was investigated in Ref. || with an arbitrary atomic 
interaction. 

For a Bose-condensed gas confined in a magnetic trap, 
the total number of particles N of the system is con- 
served and hence a canonical ensemble is appropriate to 
investigate the particle-number fluctuations in the con- 
densate. Directly from the canonical partition function 
of the system and by using a developed saddle-point 
method, a systematic app roach was proposed by the 
present authors |13l Il4l Il5l llq] for investigating the role 
of collective excitations on the particle-number fluctua- 
tions. Within the canonical ensemble a general method 
has been given recently for studying the thermodynamic 
properties of Bose-condensed gases based on the calcula- 
tion of the probability distribution function [l^. In 
Ref. EH , the theory in 0, 0] is developed to calculate 
the particle-number fluctuations due to collective excita- 
tions by including the effect of quantum depletion. 

For a 3D Bose-condensed gas confined in a mag- 
netic trap, although the effect of the quantum de- 
pletion of condensate is discussed by Giorgini et al 
6] within the particle-number-nonconserving Bogoliubov 
approach, explicit expression in the formula denoting the 
particle-number fluctuations by the quantum depletion is 
not provided. In the present work, we shall carefully cal- 
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culate the particle-number fluctuations due to collective 
excitations within a canonical ensemble. We take into 
account the effect of quantum depletion and finite-size of 
the condensate. 

We first give a brief description for the method devel- 
oped in Ref. 0] . For a condensate confined in a 3D 
magnetic trap, the collective excitations generated from 
the condensate can be described by three quantum num- 
bers n (the number of radial modes), I (the magnitude of 
the total angular momentum), and m (the z-component 
of the angular momentum). Assuming that N^ lm is the 
number of the collective excitations indexed by the quan- 
tum numbers nlm and No is the number of particles in 
the condensate, the canonical partition function of the 
system takes the form 



Z[N] =5>xp 



-0\Eo + 



/ j N nlm e n h. 



, (1) 



where the prime in the summation represents the condi- 
tion that the total number of atoms in the system should 
be conserved within the canonical ensemble, Eq is the 
energy of the condensate which can be regarded as a 
ground-state energy of the system. Using the Bogoliubov 
theory lift Il8l| and a saddle-point method developed in 
HE QUI, for the temperature below the BEC criti- 
cal temperature the normalized probability distribution 
function reads 113 
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V +Vext (r) - n + 2gn (r) + gn (r) u ni 



^nlm^nlm i 



(5) 



where V ex t (r) is the external potential confining the Bose 
gas, ji and g are respectively the chemical potential of the 
system and the coupling constant; n (r) and no (r) are the 
density distributions of the Bose gas and the condensate, 
respectively. 

For temperature far below the critical temperature, we 
have JVq >> 1. In this situation, from Eq. (J2J) the fluctu- 
ations of the condensate contributed from the collective 
excitations are given by 



(6 2 N ) = (N 2 ) - (N Q ) 



(6) 



For Eqs. an< 4 (5), u n i m and v n i m can be approximated 
as 19] 
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For a Bose gas confined in an isotropic harmonic poten- 
tial with angular frequency lo, Xnlm an( l £ nlm(= f^nlm) 
are determined by the eigen equation 



G n (N,N ) = A n exp 



2E 



(2) 



where A n is the normalization constant and Ng is the 
most probable value of the atomic number in the con- 
densate. The quantity 5 is given by 



2 V ' [(^ ~ y Xnlm] = "llrnXnlm* ( 8 ) 

where R is the radius of the condensate, determined by 
the chemical potential /i of the system through fi = 
muj 2 R 2 /2. In Eq. JfJJ, uo rdm and Xnlm are f° un d to be 
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where u n i m , v n i m and e n i m (the energy of the collec- 
tive mode nlm) are determined by the following coupled 
equations: 
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V 2 +V ex t (r) - jU + 2gn (r) u nim + gn (r) v ni 



<J->nh. 



oj (2n 2 + 2nl + 3n + l) 1/2 . 



(9) 



and 



Xnlm = A nl Pl 2n) (1) r l Y lm (0, ip) 6 {R - r) , (10) 

where A n i m is the normalization constant determined 
by J IXnimi 2 dV = 1 and (a;) is a step function. In 
Eq. (fTUfl . P ; ( 2n ) (x) = J2k=o a 2kX 2k is a polynomial with 
coefficients satisfying the recurrence relation: a2k+2 = 
-a 2k (n - ft) {21 + 2k + 2n + 3) / (k + 1) (21 + 2ft + 3) 
with ao = 1. 

With the above formulas we now calculate the particle- 
number fluctuations in the condensate. Substituting 
the above results into Eqs. J3J) and 10, we obtain the 
particle-number fluctuations due to collective excitations: 



£ nlm^nlr. 



(4) 



(5 2 N ) = »i + 9f? 2 + K 3 



(11) 
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where 



= A, (JL 



»i = Ai 



and 



3*3 
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(12) 



(13) 



(14) 



The coefficients Ai, A2, A3 are given by 



a 2 = E 



(2f+l)^ 



(2n 2 + 2nl + 3n + Z 1 ) 2 ' 



(2i + l)^ 



]~T (2n 2 + 2nl + 3n + i) 



nl=£0 



3/2 ' 



(2n 2 + 2n2 + 3n + /)' 



(15) 



with 



j? (1 - x 2 ) (p/ 2 *> (,) 



c 2l+2 dx 



Jo ( P i 2n) ( x )) x 2l + 2 dx 



(16) 



The factor (2Z + 1) in Ai, A2 and A3 is due to the 
(21 + 1) —fold degeneracy of the angular momentum. 
The numerical results of the above parameters are 
Ai = 0.47, A 2 = 0.94, A 3 = 4.96. Note that the coeffi- 
cients in the particle-number fluctuations of the conden- 
sate is different from the result given in Refs. [6j and 

0. 

In Eq. (fTTfl . 3?i is the leading term of the particle- 
number fluctuations, while K 3 shows the particle-number 
fluctuations due to the quantum depletion of the conden- 
sate which do not vanish even at zero temperature. 

For the Bose-condensed gas in the harmonic potential, 
the chemical potential is given by 



fuv /157V a^ 

2 V a hi 



2/5 



(17) 



with a being the s— wave scattering length and ato = 
yjh/mu) the harmonic oscillator length. Using the ex- 
pression of the critical temperature of an ideal Bose gas, 
T^=tkJ (N/C (S)) 1 ^ 3 /fcs, and introducing the dimension- 
less parameters t = T/T® and a — a/ahoi we obtain the 
leading term of 5Ri: 

Hi = 0.91i 2 (1 - t 3 ) 4/5 a 4 / 5 N 22 ^ 15 . (18) 



We see that the leading term of the particle-number fluc- 
tuations in the condensate exhibits an anomalous behav- 
ior of (S 2 N ) ~N 22 / 15 . 

It is straightforward to obtain the results of 3^2 and 
3?3, which are given by 



3i 2 = 1.93^1 - a 4/ V/ 5 7V 17 / 



/15 



and 



K3-2.70 (i-i 3 ) 4/ V/ 5 7v 4 / 5 . 



(19) 



(20) 



We now make some remarks on the results obtained 
above: (i)The anomalous behavior of the particle-number 
fluctuations of the condensate predicted by (fl%)) is ob- 
tained when the particle number of the system is finite. 
In the thermodynamics limit of the Bose gas in the har- 
monic potential, i. e. letting N — > 00 and to — > while 
keeping Nuj 3 constant 21], the anomalous behavior of 
the form (S 2 N ) ~ N A / 3 given in Refs. 0, El can be ob- 
tained. (ii)Different from the result obtained in Ref.0 
where only the low-temperature behavior T « Xj? is 
taken into account, in the present work the temperature 
dependence of the fluctuations is valid for the temper- 
ature region below the critical temperature. When the 
temperature T approaches the BEC critical temperature 
T c °, the particle number in the condensate approaches 
zero and hence the fluctuations originated from the col- 
lective excitations vanish, as shown in Eq. Q18J1. (iii)For a 
higher temperature the collective excitations play only a 
second role and one must consider the contribution from 
single-particle excitations. The particle- number fluctua- 
tions contributed by the single-particle excitations show 
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the behavior of (S 2 N ) - N 13]. (iv)From Eq. JTHJ , we 
see that the confining potential has the effect of increas- 
ing the particle number fluctuations. It is easy to know 
from Eq. ftjfy that 5Ri - lu 2 / 5 . 

In conclusion, we have investigated the particle- 
number fluctuations originated from collective excita- 
tions in a three-dimensional, repulsively interacting Bose- 
Einstein condensate (BEC) confined in a harmonic po- 
tential. We have carefully calculated the contribution 
to the fluctuations due to the quantum depletion and 
provided the explicit expression of the coefficient in the 
formulas denoting the particle- number fluctuations. Our 
results show that the particle-number fluctuations of the 
condensate due to the collective excitations display an 
anomalous behavior of the form < 8N§ >^ jV 22 / 15 . 
When T — » T c ° such fluctuations approach to zero and 
hence the fluctuations due to single-particle excitations 
become dominant. It is possible that the anomalous be- 
havior predicted in this work could be observed experi- 
mentally by means of, e. g., the scattering of short and 
nonresonant laser pulses on a BEC [22j |. 
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